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Abstract. We derive a simple and general relation between the fidelity of quantum 
motion, characterizing the stability of quantum dynamics with respect to arbitrary 
static perturbation of the unitary evolution propagator, and the integrated time 
auto-correlation function of the generator of perturbation. Surprisingly, this relation 
predicts the slower decay of fidelity the faster decay of correlations is. In particular, 
for non-ergodic and non-mixing dynamics, where asymptotic decay of correlations is 
absent, a qualitatively different and faster decay of fidelity is predicted on a time scale 
oc 1/5 as opposed to mixing dynamics where the fidelity is found to decay exponentially 
on a time-scale oc 1/(5^, where (5 is a strength of perturbation. A detailed discussion 
of a semi-classical regime of small effective values of Planck constant h is given where 
classical correlation functions can be used to predict quantum fidelity decay. Note 
that the correct and intuitively expected classical stability behavior is recovered in the 
classical limit ?i — > 0, as the two limits 6^0 and ?i — > do not commute. In addition 
we also discuss non-trivial dependence on the number of degrees of freedom. All the 
theoretical results are clearly demonstrated numerically on a celebrated example of a 
quantized kicked top. 
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1. Introduction 

The precise quantum signatures of different quahtative types of classical motion and 
the very definition of chaos in quantum mechanics are still the issues of an unsettled 
discussion (see e.g. ^ |^). Due to unitarity of quantum dynamics, quantum chaos 
cannot be defined in the same way as the classical chaos 0], namely through the 
exponential sensitivity on the variation of initial conditions. However, Peres |Q proposed 
an alternative concept which can be used in classical as well as in quantum mechanics: 
One can study the stability of quantum motion with respect to a small variation in the 
Hamiltonian, or more generally, a variation of the unitary evolution operator. Clearly, 
in classical mechanics this concept, when applied to individual trajectories (or to phase 
space distribution functions as we show below in sect. 4), is equivalent to the sensitivity 
to initial conditions: Integrable systems with regular orbits are stable against small 
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variation in the hamiltonian (the statement of KAM theorem), wheres for chaotic orbits 
varying the hamiltonian has a similar effect as varying the initial condition: exponential 
divergence of two orbits for two nearby chaotic hamiltonians. 

This paper is devoted to a systematic theoretical study of the stability of a unitary 
time evolution with respect to small static variations of the unitary propagator. It 
will be primarily applied to the Schrodinger propagator in quantum dynamics (sect. 2 
and 3), however an alternative application to the stability of classical unitary Perron- 
Frobenius evolution will be considered for comparison (sect. 4). The quantity of the 
central interest here is the fidelity of unitary (quantum) motion. Let us consider a 
unitary operator U being either (i) a short-time propagator U = exp{—iH At /h) , or (ii) a 
Floquet map U = T exp{—i Jq dTH{T)/h) of (periodically time- dependent) Hamiltonian 
H {H{t+p) = H{t)), ot (iii) a quantum Poincare map or any other quantized symplectic 
map. In any case a general small perturbation of the unitary operator can be written 
in the following form 

Us = Uexp{-iA6/h), (1) 

where A is some self-adjoint operator, 5 is a perturbation strength and h is an effective 
Planck constant which is taken out explicitly so that the observable A has a well 
defined classical limit (e.g. the corresponding Weyl symbol). The influence of a 
small perturbation to the unitary evolution is quantitatively described by the overlap 
{'ipsit)\ip{t)) measuring the Hilbert space distance between the exact and the perturbed 
time evolution from the same initial pure state \ipit)) = U^\ip), I'ipsit)) = Ug\ip), where 
integer t is a discrete time (in units of the period p). This defines the fidelity 

F{t) = {UiV'), (2) 

where () gives the expectation value in the initial state lip). More generally, it may 
be useful to statistically average the fidelity over an ensemble of different pure initial 
states \ipk) appearing with probabilities pk- Thus we will write the fidelity in terms of a 
statistical density operator p = J2k Pk\i^k){i^k\, namely as eq. @ using the conventional 
statistical expectation value {B) = ti pB. The theoretical discussion in this paper is 
fully general, however, we will later apply our theory in the two extreme situations, 
namely for the (coherent) pure initial state p = and for the full Hilbert space 

average p = l/M which is also equivalent to considering a random pure initial state of 
maximal information entropy. Integer M denotes a dimension of the relevant Hilbert 
space which can be written semiclassically by the Thomas-Fermi rule 

J\f = {27cn)-'^V (3) 

where V is the classical phase-space (or energy shell) volume and d is the number of 
(relevant) degrees of freedom. 

The quantity F(t), or its equivalent definitions, have already been discussed in 
several different contexts in quantum physics. We name just a few, namely those which 
provided direct motivation for the present work: First, F{t) has been proposed by Peres 
P as a measure of the stability of quantum motion. Second, |-F(t)P is the Loschmidt 
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echo characterizing the dynamical irreversibility, which has been used e.g. in spin-echo 
experiments where one is interested in the overlap between the initial state 
and an echo state U^'^U^\ip) obtained by composing forward time evolution, imperfect 
time inversion with a residual (non-invertible) interaction described by the operator 
A6, and backward time evolution. Third, the fidelity has become a standard measure 
characterizing the loss of phase coherence in the quantum computation Fourth, 
it was used to characterize the "hypersensitivity to perturbations" in related studies 
of information theoretic characterization of quantum chaos [H, though in a different 
context of a stochastically time-dependent perturbation. 

The main result of this paper is a relation of the fidelity to the ergodic properties of 
quantum dynamics, more precisely to the time autocorrelation function of the generator 
of perturbation A. Quantum dynamics of finite and bound systems has always a discrete 
spectrum since the effective Hilbert space dimension Af is finite, hence it is non- ergodic 
and non-mixing 0: time correlation functions have fluctuating tails of order ~ 1/J\f. In 
order to reach a genuinely complex quantum motion with a continuous spectrum one has 
to enforce A/" ^ oo by considering one of the following two limits: semi-classical limit 
of effective Planck's constant ^ — 0, or thermodynamic limit of number of particles, 
or number of freedoms d oo. Our result is surprising in the sense that it predicts 
the fidelity to decay slower if the integrated time correlation function is smaller, i.e. 
if the relaxation process is faster. In particular, for ergodic and sufficiently mixing 
dynamics (such that the time integrated autocorrelation function of the perturbation 
A is finite) the fidelity is found to decay exponentially on a time-scale oc whereas 
for a 'more regular', non-ergodic dynamics with non- vanishing time- averaged correlation 
functions, the fidelity decay is found to be qualitatively faster with a characteristic time- 
scale (X 6^^. However, this surprising and apparently counterintuitive result is correctly 
reconciled with the expected corresponding behavior in the classical limit due to non- 
trivial non-interchangability of the limits 6-^0 and h —* 0. In all cases, the relevant 
time-scale for the decay of fidelity in the semi-classical regime can be explicitly computed 
in terms of the classical quantities only (e.g. classical correlation functions). The main 
theoretical predictions are clearly demonstrated on the numerical examples of a kicked 
top and a pair of coupled kicked tops. Furthermore, our theory on fidelity is general 
and can be applied to any perturbed unitary evolution; as an example we consider the 
stability of the classical unitary Perron-Frobenius (Liouville) evolution for maps, where 
we show how 'classical fidelity' behaves in a qualitatively different way than the quantum 
fidelity and conforms to the corresponding linearly stable and exponentially unstable 
behaviors in the respective limiting cases of regular and chaotic dynamics. Two short 
announcements of our theory with applications in the contexts of many-body quantum 
dynamics and quantum computing have already been reported |T0, ITT . 



In section 2, we present a theoretical derivation of the general relation between 
the fidelity of unitary motion and the correlation decay, discuss qualitatively different 
regimes of stability in the semiclassical range of small effective h and outline the 
corresponding time-scales. In section 3, our theoretical predictions are demonstrated 
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in detail by numerical experiments on a quantized kicked top. In addition, we consider 
also a pair of coupled kicked tops in order to illustrate some dimensionality dependent 
aspects. In section 4, we demonstrate the conceptual difference between the stability of 
classical and quantum unitary evolution by applying our theory to the case of classical 
(Liouvillian) unitary evolution of phase space densities. In section 5, we conclude by 
pointing to some potentially important applications and implications of our results. 

2. General theory 

We start by rewriting the fidelity (0) in terms of the Heisenberg evolution of the 
perturbation At := U~^AU^. Let us recursively insert the unity f/*'?7~*' in the definition 
(0) for F{t) and observe U~*'U^^U^^^' = exp{iAt'5 /h) for t' running from t — 1 downto 
0, yielding 

F{t) = {UJ'U') = {UJ^'-^^U'-^U-^'-^\Ui^U)U'-^) = 

= {U;^'-^^U'~^exp{tAt^2S/n)exp{tAt^i6/h)) = 

= {exp{iAoS/h) exp{iAiS/h) ■ ■ ■exp{iAt^iS/h)). (4) 
The obvious next step is to expand the product into a power-series in 5 

oo -m 5:171 t—1 

^W = i + Ezni^^E (5) 

m=l ^■'^ ti,...,t^=0 

where the operator T denotes a left-to-right time ordering. Such a perturbative 
expansion converges absolutely for any 6 provided only that the perturbation A is a 
bounded operator. Therefore, the fidelity may be computed by an arbitrary truncation 
of the above expansion (^) where its order M (m < M) depends only on the desired 
accuracy (9(5^^^^). We can see that the fidelity F{t) can be expressed entirely in terms 
of time correlation functions of the generator A of perturbation. 

Before proceeding further we note that the first order term (m = 1) can be 
trivially eliminated, namely by transforming the perturbations as A A — al with 
a := (l / 1) Yft~=o{^t') (i-e. the perturbation A has to be traceless if p = l/A/"). This 
does not affect the absolute value since a shift by a multiple of unity reflects 

in a simple complex rotation of fidelity F{t) F{t) exp{—i6at/h). So, to the lowest 
non-trivial order (M = 2) we find 

F{t) = l-^ E fC{t'X) + 0{6% (6) 

t',t"=0 

where 

C{t', t") := {At'Atn), fC{t', t") = C(min{t', t"}, max{t', t"}), (7) 

is a 2-point time correlation function of the quantum observable A. In the following, the 
regime of sufficiently small 6, such that the fidelity remains close to 1 (i.e. 1 — F{t) <C 1) 
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and the above formula is accurate, will be referred to as the regime of linear response. 
Although we have a general convergent formula for systematic higher order corrections 
(|) we note that all the qualitative and in most cases even quantitative results of this 
paper can be explained by this simple linear- response Kubo-like formula (^). 

An interesting and somewhat counterintuitive conclusion can be drawn from the 
expression (^), namely the smaller the time-integrated correlation function the higher 
the fidelity. In the semiclassical regime of approaching the classical limit h ^ the 
quantum correlation function C{t', t") goes over to the corresponding classical correlation 
function provided a state p converges to some classical phase space distribution, like e.g. 
for p = 1/U. 

Therefore, the fidelity for classically chaotic systems will decay with the rate which 
is inversely proportional to their rate of mixing, and furthermore for classically non- 
ergodic, i.e. regular or integrable motion, the correlation functions will generally not 
decay to zero and the fidelity will therefore decay much faster. Of course, this surprising 
general conclusions can only be true if the initial state (or the statistical density 
matrix p) is sufficiently random., i.e. if it is not a small linear combination of the 
eigenstates of U (or Us) as will be discussed in more detail later. Now let us have a 
closer look at fidelity in the two qualitatively different regimes of correlation decay. 

2.1. Regime of ergodicity and fast mixing 

Here we assume that the system is (classically) ergodic and mixing such that the 
correlation function of the perturbation A decays sufficiently fast; this typically (but 
not neccessarily!) corresponds to globally chaotic classical motion. Due to ergodicity 
we will safely assume p = i.e. average over the full Hilbert spacefj] (.) = tr(.)/A/'. 

For any other initial state (e.g. in the worst case for the minimal wavepacket - coherent 
state) one obtains identical results on F{t) for sufficiently long times, i.e. longer than the 
Ehrenfest - ergodic time ~ ln(l/^)/A (for a classically chaotic system with maximal 
Lyapunov exponent A) needed for a minimal wavepacket to spread effectively over the 
accessible phase space. The state averaged quantum correlation function is homogeneous 
in time C{t\ t") = C{t" — t') so we simplify the second order - linear response - formula 
for the fidelity 

m = 1 - ^ |itc(o) + 1: {t - t')c(t')} + o{6'), (8) 

and note again that we have assumed A to be traceless (A) = 0. 

If the decay of correlation function C{t) is sufficiently fast, namely it is sufficient 
that the sums Ik = St^-oo exist for k G {0, 1} and that a certain characteristic 

mixing time exists, e.g. tmix = h/^o, then the above formula can be further simplified. 
For times t ^ t^i^ we can neglect the second term under the summation in (^) and 

II In case of autonomous (time-independent) systems, canonical or micro-canonial state should be used 
instead. 
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obtain a linear decay in time t in the linear response regime 

F^t) = 1 - (S/nyat with a = -C(0) + ^ (9) 

^ t=i 

which is a transport coefficient given by the integral of correlation function. We can 
make a stronger statement in a non-linear-response regime if we make an additional 
assumption on the factorization of higher order time-correlations - n— point mixing. 
This implies that 2m-point correlation (A^^ ■ ■ ■ At^^) is appreciably different from zero 
for t2m — ti — i> oo only if all the (ordered) time indices {tj,j = 1 . . . 2m} are paired with 
the time differences within each pair t2j — t2j-i being of the order or less than tmix. Then 
we can make a further reduction, namely if t ^ rnt^^^ 

r Y {KK---KJ^T Y (^*.A.)---(A,_A,j^i-i:(ta)- (10) 

tl,---,t2,n=0 ti,...,t2m=0 

SO that we obtain a global exponential decay 

Fera{t) = exp {-t/r^^), T^^ = — , (11) 

with a time-scale rem = 0{6^^). We should stress again that the above result (|TI|) has 
been derived by the assumption of true quantum mixing which can be justified only in 
the limit JV —>■ oo, e.g. either in semiclassical or thermodynamic limt. However, precise 
study of finite size effects and the relevant time and perturbation scales will be given in 
subsection 0|. 



2.2. Non-mixing and non-ergodic regime 

The opposite situation of a non-mixing and non-ergodic quantum dynamics, which 
typically corresponds to integrable, quasi-integrable (KAM), or mixed classical 
dynamics, is characterized by a non-vanishing time-average of the correlation function 

C= limi Y fC{t\t"). (12) 

t',t"=0 

Here, due to non-ergodicity, time-average C depends on the choice of the initial state 
or, more generally, on the density matrix p. We assume that a characteristic 
averaging time-scale tave exists, namely it is an effective time t = tave at which the 
limiting process ([T^) converges. Therefore, for sufficiently large times t ^ tave the 
linear-response formula @ gives, in contrast to (P), a simple quadratic decay in time 

1 / t \^ h 



F„e(t) = l-^f— ) +0{S'), rne = -^, (13) 

with time-scale Tne = 0{6~^). One should observe that the non-ergodic time-scale Tne 
can be much smaller than the ergodic-mixing time-scale Tem ([Tl|) provided h is fixed, or 
the limit 5 — > is taken prior to the limit ^ — >■ 0. Again, we can make much stronger 
general statement using an additional assumption. Namely, if we assume that t ^ tave 
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then we can re-write the m-tuple sums in the series in terms of a time average 
perturbation operator 

A = ^lim(l/t)|] A,, (14) 

t'=0 

namely 

Note that A is by construction an integral of motion |T^, [U, A] = 0, and reduces to a 
trivial multiple of identity (in fact, it vanishes with our choice of (A) = 0) in the ergodic 
case of previous subsection. Whereas in an ergodic and mixing case, m— th order term 
of (^ grows with time only as (9(t™'/^) (for even m) since it is dominated by pair time 
correlations, here in a non-ergodic case, the non-trivial time average operator A already 
gives the dominant effect, namely 0{t"^) for m— th order term of (^, so the effect of pair 
time correlations can safely be neglected for sufficiently long times (t ^ tave)- Observe 
also that time averaged correlation is just a second moment of the time average operator 

C={A'). (16) 

In the semiclassical regime, C goes to a purely classical, /i-independent quantity 
only for such initial states {ip), corresponding to /i-independent distribution of initial 
conditions in the classical phase space, e.g. for random states or p = Formula (^) 

expresses the fidelity in terms of a sequence of moments {A"^). It has been shown [1^ 



that for a general class of one-dimensional spin chains, the moments in thermodynamic 
limit d oo oi long chains tend to a normal behavior (A^™) (2m — 1)!!(7™', 
i/^2m+i\j _^ fidelity exhibits a gaussian decay 

i^normal(t) = exp {- {t / T^,f / 2) . (17) 

Although good numerical agreement with this formula has been observed in several 



related contexts, e.g. in quantum computing |TT], it cannot be generally valid, in 
particular not in the semiclassical regime with a small number of freedoms d, as we 
show below. 

In general, we can rewrite the formula (|T3|) in terms of (average) local density of 
states (LDOS) of the time average A which is defined (for finite M) as 

dp{a) = ^5{a- an){an\p\an). (18) 

n 

where a„ are the eigenvalues and |a„) the eigenstates of self-adjont operator A. Namely 

{e^^(i5At/h)) = Jdae'^'^/^dpia), (19) 

so the fidelity -Fne(^) of non-ergodic systems is just a Fourier transformation of LDOS 
dp{a) of the time average perturbation A (which becomes the simple (global) density of 
states in case of an uniform average p = l/M). Thus only a gaussian LDOS would yield 



a gaussian fidelity decay (17), however we often find good finite time gaussian behavior 



Stability of quantum motion and correlation decay 



8 



which goes beyond the second moment Jdaa'^dp{a) = (A^) = C, even if the asymptotic 
time dependence is completely different, see e.g. fig. ^for an example of kicked top. 

Let us now try to give some semi-classical estimates on the fidelity decay in non- 
ergodic regime by expressing the asymptotics ^ — >^ in terms of classical information 
only. Here we assume that our (classical) map is integrable so that there exist d canonical 
constants of motion - the action variables / = (Ji, . . . , J^), which are quantized using 
EBK rule /„ = h{n + 7/4) where integers Uj > 0, j G {1, . . . , rf} are quantum numbers 
and integers 7j Maslov indices. Since the time average operator A commutes with U 
and with the actions /, it is diagonal in (generically non-degenerate) basis of eigenstates 
of I (quantized tori) |n), namely in the leading semiclassical order 

{n\A\n') = 5^,r^'d{Ir^) (20) 

where a(/) is the corresponding classical time-averaged observable in action space. The 
fidelity (|l^) can therefore be written as 

Fne{t) = ^exp(zta(/„)V^)(n|p|n) (21) 

n 

Now, provided the diagonal elements of the density matrix can be written in terms of 
some smooth function Dp{In) = {n\p\n), and replacing the sum (^T]) by an integral 
over the action space, which is justified for small h up to classically long time oc h^5~^, 
we obtain 

F^,{t) = h-'^ J d'^I exp {ita{I)6/h}Dp{I) (22) 

The obvious next step is to compute this integral by a method of stationary phase. 
However, the result depends on the precise form of the function Dp{I) (which may also 
explicitly depend on h) so we only work out the details for two important special cases. 



2.2.1. Semiclassical asymptotics for an uniform initial state average. Let us first 
assume uniform averaging over (random) initial states p = l/A/", so Dp{I) = 1/M = 
(27ihY/V. Then, for large t6/h the above integral ( p2D can be written as a sum of 
contributions stemming from, say p points, J,,, ?7 = 1, . . . , p where the phase is stationary, 
dd{ljj)/dl = 0, yielding 



Frit) = 



'^^^ V" exp{ita{Q5/h + iu^} 

^ 'detAJV2 ' ^^^> 



where {A^j^fc := d'^d{Ir,) / dljdlk is a matrix of second derivatives at the stationary point 
7], and = 7r(m+ — m_)/4 where m± are the numbers of positive/negative eigenvalues 
of the matrix A^. We repeat that the stationary phase formula (|23| ) is expected to be 
correct in the range const h/6 < t < const'/(5. Most interesting to note is the asymptotic 
power-law time and perturbation dependence F^^"^ ~ \h/{t5)\^^'^, which allows for a 
possible crossover to a gaussian decay (p!7D when approaching the thermodynamic limit 
d ^ oo. 
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2.2.2. Semiclassical asymptotics for a coherent initial state. Now let us consider a 
single (i-dimensional general coherent state centered at {I*,0*) in action-angle space, 

p = \i*,e*){r,e*\, 

{n\r, = f - j |det A|^/Sxp |-^(/r^ - /*) ■ A(/^ - r) - m • ^ | , (24) 

where A is a positive symmetric d x d matrix of squeezing parameters, giving Dp{I) = 
{h/7rY/^ \det A\^^^ exp{- {I - r) ■ A{I - r)/h) and 

Idet M^^^ r r 1 1 

F-'^(t) = J d'^i exp |- {^t-a{I)S -{I- I*) ■ A{I - I*))] ■ (25) 

Using the assumption 6t <^ 1, we see that a unique stationary point Is of the exponent 
approaches /* as 5 0, 

I^ = r-'-^A-'a' + Oi5^), where a' := (26) 



so we may explicitly evaluate ( pS]) by the method of stationary phase without any lower 
bound on the range of time t, 

^o.W.exp{-(i^^4^.^.!^| (27) 

Note that the fidelity decay for a coherent initial state with regular classical motion has 
a time-scale 

Tne-coh = (2/i)'/'(a'- A-^a')-'/'(5-^ oc (28) 

which is consistent with ( pTSD with C = ^h{a' -A^^a') and is typically (for small h) larger 
than the time-scale T^e c<. h/6 for a random initial state. However, surprisingly enough, 
for sufficiently small 6, fidelity decay for coherent state with regular classical dynamics 
can still be faster than the fidelity decay for ergodic and mixing dynamics. It should 
be noted that the above derivation for fidelity decay of a coherent state ( p^ - |28| ) remains 
valid in a quasi-integrable (KAM) situation of mixed classical phase space, provided 
that the initial wave packet is launched in a regular region of phase space where (local) 
action-angle variables exist. 

From the derivation of gaussian behavior (^) we can also begin to understand the 
reason for a faster gaussian decay of fidelity for a non-mixing system in contrast to the 
exponential decay for a mixing case. Namely, in the non-mixing case, the displacement 
between centers of the wave packets of two nearby regular evolutions grows hallistically 
and thus produces destructive interference (i.e. decay of fidelity) much faster than in 
the case of mixing dynamics where the destructive interference builds up in a diffusive, 
random way. 

2. 3. Finite size effects and time and perturbation scales 

The theoretical relations of the previous subsections are strictly justified in the 
asymptotics M ^ oo. However, in any experimentally or practically relevant situation 



Stability of quantum motion and correlation decay 



10 



the dimensionahty of the relevant Hilbert (sub)space is finite, so we have to discuss 
in detail the limitations and additional time and perturbation scales associated with 
this. 

First of all, for finite Af, fidelity F(t) cannot decay indefinitely but starts to fluctuate 
for long times due to discreteness of the spectrum of the evolution operator U. Let us 
write the eigenphases of U and Us, and the corresponding eigenvectors, respectively, as 
0„, 0^, and 10^), n = l,...,Af, satisfying f/|0„) = e-^-^-l^,), f/,|0^) = e-^<^"|0^). 
Let us define a unitary operator V which maps the eigenbasis of U to the eigenbasis 
of Us, namely V\(j)n) '■= \(pn) ^ill n, with matrix elements Vmn '■= {(pm\V\<Pn) , and 
write the matrix elements of the statistical averaging operator as pmn '■= {(pm\p\4>n) ■ 
Note that the matrix Vmn is real orthogonal, V G 0{Af), if U and Us possess a common 
anti- unitary symmetry (e.g. time-reversal). Now it is straightforward to rewrite the 
fidelity (|) as 

Fit) = Y.ipV)niV:i exp - 0f )t) , (29) 

n,l 

where {pV)ni = Y^mRnmYmi IS a matrix element of the operator product pV . At this 
point we are interested in the long time fluctuations so we compute the time average 
fidelity fluctuation 

FLs--=W =T.\ipy)nl?\Vnl?. (30) 

n,l 

In the process of averaging over the time we have assumed that the eigenphases are 
non-degenerate so exp(i(0„ — 0„/ + 0f — 0f')i) = 5nn'5i,i'- So we see that fidelity fluc- 
tuation Fj-ms depends on the orthogonal/unitary matrix Vmn and the statistical matrix 
Pmn- Clearly, Vmn converges to an identity matrix for sufficiently small 5, while for large 
enough 5 and sufficiently complex perturbation A one may assume that it becomes a ran- 
dom orthogonal/ unitary matrix. Of course, the critical perturbation strength 5 = 5rm 
for a transition between the two regimes is very much system dependent and cannot be 
discussed in general. Let us now discuss three different limiting cases of the statistical 
matrix pmn- 

(i) First, consider the simplest case where the initial state is a pure eigenstate of U 
say, p = |0i)(0i|. Then pmn = ^m.i^n,! and ( pO]) rewrites to 

FLs-: = j:\yu\'. (31) 

Now, for weak perturbation 6 <^ 5rm fidelity does not decay at all -Frnis'-i — 1' while 
for strong perturbation S ^ 6rm we may estimate (pID using the fact that Vmn become 
real/complex gaussian random variables with variance l/A/" for the orthogonal(/5 = 
l)/unitary(/3 = 2) case, namely FrmT-f = [(4 — P)/-^]^^"^- This result can be general- 
ized to the case where initial state can be written as a finite superposition of a small 
number, say K, of eigenstates of U: then for a weak perturbation F^f^ ~ K~^^'^. With 
this simple observation we can easily explain the numerical result of Peres [H where 
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no-decay of fidelity was found for a coherent initial state sitting in the center of elliptic 
island, thus being a superposition of a very small number of eigenstates (it is almost an 
eigenstate), while the behavior in generic cases may be drastically different as described 
in the present work. 

(ii) Second, consider the case of a random pure initial state \ip) = l]„c„|0„), giving 
Pmn = CmC* . The principle of maximal entropy forces the coefficients c„ to be indepen- 
dent random complex gaussian variables with variance l/N'. Thus, averaging (^) over 
c„'s yields 

FLs-, = ^ \Vnl\' + E \Vnl\'\V^l\'] . (32) 



n,l n,m,l 



Here, for weak perturbation 5 <^ 5rm we find -Frnil^ii = {2/J\fy/'^ while for strong per- 
turbation 6 3> (5rm we have Fj-mT-fi = M^^^'^ 

(iii) Third, in the limiting case of a uniform average over all initial states p = l/A/", 

Pmn = Smn/^f, WC have 

-^mis-iii = T7^ E l^"'!"^- (33) 
■'^ n,l 

Again, for weak perturbation 5 <^ 5rm, -^ms'-iii = A/"^^/^, and for strong perturbation 

'5»5rm,F;r-ii = (4-/5)VVAr. 

Note that the formulae (pT]-^) state simply that fidelity fiuctuation is an inverse par- 
ticipation ratio (IPR) of the perturbed eigenstate(s) (single one in case (i) or an average 
over all in case (iii)) in terms of the unperturbed eigenstates, and is thus directly related 
to the localization properties of eigenstates of Us in terms of eigenstates of U. However, 
except for the pathological case of initial state being a small combination of eigenstates 
of U with weak perturbation, the fidelity fiuctuation is always between the limiting 
values A/"^^/^, and ^/3/J\f. For an illustrations of transitions between weak and strong 
perturbations in case (iii) of a kicked top, see fig. ^ 

Therefore, fidelity will decay according to one of the asymptotic formulae 
(0,0,^,0) until it reaches the value of finite size fiuctuation. This condition 
determines the time-scale = t*(A/'), 

Fv-_oo(t*(Ar))=Fr^,(Ar). (34) 

In non-ergodic and non-mixing regime this is the only relevant finite size time-scale, 
while in the mixing regime the situation is more complicated. 



2.3.1. Regime of ergodicity and fast mixing. Here, combining (|3^) with an exponential 
decay ([TTD we get 

t, = pr,^lnAf^^Hl/h) (35) 
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Figure 1. Finite-size fidelity fluctuation for a kicked top ( pSf ) in the mixing regime 
(a — 30) and with the whole Hilbert space average ( p3| ) is shown as a function of the 
perturbation strength S at four different values of J(= l/K) — 100, 200, 400, 800 (curves 
from right to left, respectively). Data are scaled in such a way J^^^-Frms (note that 
M = J) that the weak perturbation regime corresponds to value 1 while the horizontal 
lines give the corresponding values in the strong perturbation regime F'rms'-fii ~ V^/ 
Note that the crossover perturbation scales as « 1/J. 



where /z = — In Frms/ InA/" is a parameter which typically lies in the range 1/2 < < 1 
(except if the initial state is non-random and the dynamics is non-ergodic (e.g. case (i)) 
and the perturbation is small 5 <^ 5rm) and depends on IPR and the statistical operator 
p as discussed above, and a^i is the classical limit of the transport coefficient (^. 

The second new time-scale is related the to asymptotic non-decay of time 
correlations for finite M quantum dynamics, namely even if the system is classically 
mixing the quantum correlation function will have a small non-vanishing (^-dependent) 
time average 



C={A{t)A{t'))=Y.PnnAl^ (36) 

n 

where Amn = (0m|^|0n)- However, since the classical system is ergodic and mixing, we 
will use a version of quantum chaos conjecture saying that ||T^ then A^n are independent 
gaussian random variables with a variance given by the Fourier transformation S{uj) of 
the corresponding classical correlation function Cc\(t) at frequency uj = 0m — 4>n- On 
the diagonal we have a; = and additional factor of 2 due to random matrix measure 
on the diagonal (see e.g. IH), so we can write 

C=^^, (37) 
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where = S{0)/2 is the classical limit of (^. Due to ergodicity, for large A^, this 
does not depend on the statistical operator p; neither in cases (ii,iii), nor in case (i) if 
one assumes additional ensemble averaging. The decay of fidelity will start to be 
dominated by the average plateau at time t, when J2l~2o {t - t')C ^ laJ^jM > a^t, 
i.e. for times t greater than tp 

tp = IaT oc h-"^ (38) 

which is just the Heisenberg time associated to the inverse density of states. 

Now, depending on the interrelation among four (or five) time-scales rem oc h'^5~'^, 
tp oc h-'^S°, t, oc hHn{l/h)S-'^d, tmix oc h^'^, (and oc \n{l/h)S'^ if we are considering 



coherent initial states, like e.g. [^4|, [T^, |T6[) we can have four (or five) different regimes 
depending on three main scaling parameters: perturbation strength 5, Planck's constant 
h, and dimensionality d. Note that we always have t^, > rem- All different regimes can 
be reached by changing only 6 while keeping h and d fixed (see fig. 0): 

(a) For sufficiently small perturbation 6 we will have tp < rem. This means that 
-^em(^p) ~ 1 and we will have initially quadratic decay (|13]) with C given by an av- 
erage finite size plateau (^7\). This will occur for 6 < 6p where 



In fact, in this regime, also referred to [0, as perturbative, one may use a first order 
stationary perturbation theory on the eigenstates of Us, yielding 0^ = 0„ + Ann5/h + 
0{S'^), Vmn = Smn + 0{6), and rewrite (following ^7^) the finite size fidelity ( pOj) in terms 
of a Fourier transform of a probability distribution w{A) of diagonal matrix elements 
Ann, Fp(t) = J dAw{A) exp{—iAt6/h). Since w{A) is conjectured to be gaussian for 
classically ergodic and mixing system [jl3[, it follows that Fp(t) is also a gaussian with 
a semiclassically long time-scale Tp 

f,M=exp(-(<A,)V2). .,= (-) - = ^^-^— . (40) 

(b) If < tp < t* we will have a crossover from initial exponential decay of fidelity 
(P!T|) to a gaussian decay (PDj ) at t ~ tp, which will terminate and go over to fluctuating 
behavior when -Fp(t) = Fi.ms(A/'). Note that this will happen before time t^, which is 
estimated (|^) based on a slower exponential decay (0). This regime will exist in 
perturbation range 6p < 6 < 6s with an upper border 6s determined by the condition 
tp t^ 

4 = (plnAr)i/\. (41) 

(c) If we still further increase 6, we have the most interesting, 'fully nonperturbative' 
regime, when rem < t,,, < tp and we will have a full exponential decay ([TT|) , up to time 
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Figure 2. Schematic view of different regimes (a-d) of fidelity decay. 



(42) 



when the fidehty reaches finite size fiuctuations. This regime continues for 5 < 5n 
where the border 

"mix — r—; I „ I Op- 

V cl''mix ^^'■mix^ 

is determined by the condition = t^bi which is a point where the arguments leading 
to the factorization (^) and exponential decay ( |TT|) are no long er valid. We note that 
the relative size of this window range 5mix/5s = \Jj^/ (^mix In A/") increases both, in the 
semiclassical and in the thermodynamic limit. This regime also corresponds to 'Fermi 
golden rule decay' discussed in |16[ . 



(d) Further increasing 5 > S^m, the estimated fidelity decay time eventually becomes 
smaller than the classical mixing time tmix, Tem < ^mix- In this regime, the perturbation 
is simply too strong so that the fidelity effectively decays within the shortest observable 
time-scale (tmix)- 

However, if we consider a non-random, e.g. coherent initial state then the quantum 
correlation function relaxes on a slightly longer, namely Eherenfest time-scale so the 
regime (c) should terminate already at a little smaller upper border 6 = 6-^ which is 
naturally determined by the condition r^m = 

^ [adln(l//i)]i/2 [ln(l/;i)]V2- ^ ' 

For coherent initial states one thus obtains an extra but very narrow regime 5e < S < 
fornix (describing the time-range tmix < t < t-^) where the fidehty decay can be computed 
in terms of classical Lyapunov exponents |ll4| , |15|. 

2.3.2. Non-mixing and non-ergodic regime. In this regime things are simpler, as we do 
not have to worry about the average plateau in the correlation function due to a finite 
N because we already have a higher average time correlation C Cd 7^ ( [T^ ) due 
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to non-mixing nature of dynamics. Thus we have here only two relevant time-scales, 
namely Tne giving initial quadratic decay ([T3|), and the finite size fiuctuation time-scale 
t* which depends on the properties of the initial state (power law ( p3D for a random 
initial state, versus gaussian ( p?]) for a coherent initial state). We conjecture, based on a 
rigorous result in spin 1/2 chains ||10|, that the fidelity decay in the thermodynamic limit 
d ^ oo will generically approach a gaussian (|T^, which is consistent with increasing 
power d/2 oi power- law fidelity decay (|^) for finite dimension d. 



We can summarize our findings by the following question: Can we have a regime with 
Tjic < Tern, SO that the fidelity will decay faster for a "regular" (non-ergodic) than for a 
"chaotic" (ergodic and mixing) dynamics? Indeed, if 

^ ^ f Sj- := hClC^ /ac\ oc (5mix oc h random initial state 

\ 5c:=h^/'^{a' ■ K-^a'f/'^ /{^/2a^i) (X h^'"^ coherent initial state 

where ad and Cc\ are computed for a mixing and non-ergodic classical dynamics, 
respectively, then the decay in a non-ergodic case will be faster than in a mixing case. 
We can see that the condition can be generally satisfied for a random initial state, 
while for a coherent initial state it can be satisfied above the (second) perturbative 
border (see fig. |^), 5c > 6 > 6s, only in more than one dimension d > 1. We note 
that our result is not contradicting any of the existing findings on quantum-classical 
correspondence. For example, a growth of quantum dynamical entropies [18| persists 
only up to logarithmically short Ehrenfest time t-^, which is the upper bound for the 
validity of another semiclassical approach to the fidelity decay 0, jTSf and within which 
one would always find F^f'it) > Fcmit) above the perturbative border 6 > 6p, whereas 
our theory reveals new nontrivial quantum phenomena with a semiclassical prediction 
(but not correspondence!) much beyond that time. If we let ^ ^ first, and then 
6^0, i.e. we keep S ^ 5r,c(^)5 then we recover a result supported by a classical 
intuition, namely that the regular (non-ergodic) dynamics is more stable than the chaotic 
(ergodic and mixing) dynamics. On the other hand, if we let 5 ^ first, and only after 
that h ^ 0, i.e. satisfying (H), we find somewhat counterintuitive results saying that 
chaotic (mixing) dynamics is more stable than the regular one. We can conclude the 
section by saying that we have three non-commuting limits, namely the semiclassical 
limit h the perturbation strength limit 6^0, and the thermodynamic limit d —>■ oo, 
such that no pair of these limits commutes. 



3. Numerical experiments 
3.1. Kicked top 

Here we wish to verify and demonstrate the results of the previous section by numerical 



experiments. For this purpose we choose the Haake's quantized kicked top |jT9[ since 



this model served as a model example for many related studies P, |T6|, |1^, ^ |22|, [23[|. 
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The unitary propagator reads 

U = U{a, 7) = exp {-i'yJy) exp {-iaJ^/2J), (45) 

where Jk {k = x, y, z) are quantum angular momentum operators, [Jk,Ji] = i^kirJr- 
(Half)integer J determines the size of the Hilbert space 2J + 1 and the value of effective 
Planck's constant h = 1/J. The perturbation is defined by perturbing the parameter 
a, Us = U{a + 6, 7), so that the generator A is 

^ = U4V. (46) 



2 V J 

Physically, the system ( p5D represents a twist around z-axis followed by a rotation for 
an angle 7 around y-axis. The classical limit is obtained by letting J = 1/h 00 
and writing the classical angular momentum in terms of a unit vector on a sphere 
r = {x,y,z) = J /J. The Heisenberg equation for the SU(2) operators J, J' = WJU, 
reduces to the classical area preserving map of a sphere 

x' = cos 7(0; cos az — y sin az) + z sin 7 
y' = y cos az + X sin az 

z = 2; cos 7 + sin 7(2/ sin — x cos az). (47) 
Note that in the classical limit the perturbation generator is 

a{r) = ^. (48) 

For a = the system is integrable, while with increasing a there is a transition to 
chaotic motion. The second parameter 7 is usually set to 7r/2, however in our numerical 
simulation we will use two different values exhibiting qualitatively different correlation 
decay (for large a): the 'standard' case 7 = 7r/2 where Cdit) decays in oscillatory way 
and the case 7 = 7r/6 where Cdit) decays monotonically (see fig. ^. 

In the case of 7 = 7r/2 we have two discrete symmetries. The evolution U commutes 
with i?x and Ry, the rotations for vr around x and y axes, respectively. The Hilbert space 
is therefore reducible into three invariant subspaces (using notation of Peres's book 
with the basis \m) of eigenstates of and assuming J to be an even integer): EE of 
dimension J/2 + 1 with the basis states |0) and {|2m) + | — 2m)}/-\/2; 00 of dimension 
J/2 with the basis {|2m — 1) — | — (2m — l))}/-\/2; OE of dimension J with the basis 
{|2m)-|-2m)}/v/2 and {|2m- 1) + 1 - (2m-l))}/y2 with m in all three cases running 
through m = 1, . . . , J/2. For 7 7^ tt/2 the spaces 00 and EE coalesce as Ry is the only 
discrete symmetry left. In numerical experiments we always choose OE subspace so that 
the dimension of the Hilbert space is Af = J. 

We will compute the fidelity using three different statistical operators p: (1) 
p = 1/Af corresponding to full Hilbert space average, (2) pure random initial state 
P = (components Cm = {m\ip) being independent gaussian pseudo-random 

numbers) giving the same results as (1) however with a higher finite size fluctuation 
plateau Frmt (as discussed in subsect. ^73|) , and (3) pure minimal wavepacket initial 
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state p = \'d,ip){'d,ip\, namely SU(2) coherent state \^,ip) centered at the point 
n = (sin cos sin i) sin if, cos i)) on a unit sphere 

l^'¥') = E (x ) cos (79/2)-^+™ sin (^?/2)^-™e-''"^|m). (49) 
m=—J ~\~ ni / 

3.1.1. Mixing regime. In this regime we will demonstrate two main different decays 
of fidelity as discussed in subsect. pl3| , namely the exponential decay Fcmit) ([TT| ) for 
5s < 5 < (5niix and a finite size (perturbative) gaussian decay F-p{t) (^0]) for a small 
perturbation 5 < 5^. We choose large a = 30 to ensure fast mixing. This value seems to 
be excessively large, but for smaller a we still have "sticky" structures in classical phase 
space which in turn cause slow algebraic tails in the classical correlation function. We 
must remember that the derivation of an exponential decay required well defined mixing 
time-scale tmix which is not the case if we have a slow power-law decay of correlations. 

Exponential decay time ( P^ is determined by an integral/sum of time 
correlation function (^), which can be calculated in the semiclassical regime (ad) by 
means of the correlation function of the classical map 

C^,{t) = \{z\^yz\t))- (50) 

= — (z^) is a "traceless" perturbation where averaging over the sphere gives 
(z^) = 1/3. The classical correlation functions calculated by using the classical map 



(47) are shown in fig. ^. For 7 = 7r/2 the correlation function is oscillating with an 
exponential envelope hence the transport coefficient ad = 0.00385 is quite small. For 
7 = 7r/6 the correlation decay is monotonic and exponential with ad = 0.0515. The 



decay of quantum fidelity (11) can now be obtained by using the c/ass«ca/ limit a — > ad- 
Fe^(t) = exp(-5Wt), (51) 

This formula has been compared with the exact numerical calculation of fidelity (0) 
where averaging over the whole Hilbert space has been employed, i.e. p = 1/J\f, and, 
as expected due to ergodicity, there was no difference for sufficiently large J when we 
have chosen a fixed coherent initial state. As the finite size fidelity fluctuation level 
Frms-iii ( PB| ) decreases with increasing Hilbert space dimension, we chose large J = 4000 
in order to be able to check exponential decay (|5lD over as many orders of magnitude as 
possible. The results are shown in fig. |. The smallest and the largest 6 shown roughly 
correspond to borders 4 and S^i^, respectively. As we can see, the agreement with an 
exponential decay is excellent, at least over four decades in probability |F(t)p. Note 
that for 7 = 7r/2 and the largest 6 = 3- 10^^ the time-scale of the decay of fidelity 
is comparable to the classical time-scale t^i^ so the factorization assumption ([T0|) is 
strictly not applicable any more. However, due to oscillatory correlation decay, overall 
agreement with the theory (|5TD is still rather good but the oscillations in the correlation 
decay reflect in the oscillations of the fidelity decay (around the theoretical exponential 
curve). Of course, one does not need such a large J in order to have an exponential 
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Figure 3. Classical correlation function Cc\{t) ( |50| ) for a — 30, and 7 = 7r/6 (top, 
full curve) and "f ~ tt/2 (bottom, broken curve). Finite time integrated correlation 
function is shown in the inset, converging to CTcI — 0.00385 and 0.0515, for 7 — 7r/2 
and 7r/6, respectively. Averaging over 10'^ random initial conditions on a sphere is 
performed. 



decay, but for smaller J the fluctuation level Fy-y^s will be higher so the exponential decay 
(0) will persist for correspondingly smaller time, namely up to )f:*(J) (|3^ . 

Then we focus on the so-called perturbative regime S < 6p where the fidelity decay 
will be dictated by a finite size correlation average (p7|), so according to eq. ( ^OD 

Fp{t) = exp (-25V(Tcit^). (52) 

We numerically computed C (|12]) for J = 1500, a = 30 in order to show that it is given by 
the theoretical value (|3^). Quantum correlation function has been computed C{t',t") = 
{AfAt") by means of a traceless perturbation A = \{Jz/ JY — ^[(2J+ 1)(J+ 1)/J^]1. 
In fig. ^ we show a finite time correlation sum a{t) = ^ Ylt^t"=o C(t', t") which exhibits 
a crossover, at the Heisenberg time tp = J/2, from the plateau given by ad to a 
linear increase Ct due to finite size correlation average (|37D C = Aad/J- The excellent 
agreement between prediction ( ^2] ) and full numerical calculation of fidelity is shown in 
fig. 1^. In view of our findings this so-called |T^ perturbative regime can be understood 
as a simple consequence of a finite Hilbert space dimension. For times larger than the 
Heisenberg time tp every quantum system behaves effectively as an integrable one, e.g. 
with a finite time average correlation plateau. 

3.1.2. Non-mixing regime. In a kicked top this regime is realized for sufficiently 
small value of a. If the classical phase space has a mixed (KAM) structure, non- 
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Figure 4. Quantum fidelity for a kicked top with parameters a = 30, J = 4000 and 
full trace average. Top figure is for 7 = tt/2 and for 5 = 5- lO^"', 1 • 10^'^, 3 • 10"'^ 
(solid, dashed, dotted curves, respectively) Bottom figure is for -f — tt/G and S = 
1 • 10~^,2 • 10~^,3 • 10~^, (solid, dashed, dotted curves, respectively). Chain line in 
both cases gives the theory (^ij) with classically computed ad- Note that the largest 
6 — 3 ■ 10~^ case in the top figure (dotted curve) corresponds to r w 2, so it is already 
over the upper border of the regime (c, subsect. 2.3) 5 > (5mix but the agreement 
with the theory ( |5l| ) is still quite good, appart from oscillations. This is due to 
the oscillatory nature of time-correlations making the factorization assumption ( p^ 
justified (on average) even for much smaller time t as required. 
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Figure 5. Finite time quantum correlation sum G{t) = Et'Tt"=oC(*''*")/2i (solid 
curves) and the corresponding classical sum CTci(i) — Yl\^t"=o'-'c\{t' ,t")/'^t (dashed 
curves saturating at CTcI and ending at t ^ 1000) for a = 30, J = 1500. Upper curves 
are for 7 — 7r/6 while lower curves are for 7 — 7r/2. Chain lines are best fitting 
asymptotic linear functions corresponding to Ct, Q.0077t/J for 7 = 7r/2 and 0.103t/J 
for 7 — 7r/6. 

mixing regime of fidelity decay may be obtained by choosing a localized initial state 
(e.g. coherent state) located in a regular part of the phase space. However, such a 
situation may easily lead to the opposite conclusion (as compared to generic situation) 
for insufficiently large dimension Af. As discussed in subsect.^^, the fidelity fiuctuation 
plateau is determined by the number of constituent propagator eigenstates which 
are effectively needed to expand the initial state. For a coherent state sitting inside a 
(not too large) regular (KAM) island this number can be fairly small for numerically 
realizable Hilbert space dimensions, thus prohibiting any significant fidelity decay as 
observed in p. Nevertheless, we would still see the initial quadratic decay in the linear 
response regime but we would not be able to verify higher orders in the long-time 
expansion of fidelity. Therefore, in order to make a situation numerically as clean as 
possible, we choose a small value of parameter a = 0.1, such that the classical dynamics 
is almost integrable and that the majority of phase space corresponds to regular motion 
so that the number of constituent eigenstates for coherent states is as large as possible 
(on average). 

Here we focus only on the case 7 = 7r/2. For small a, the quantum and classical 
evolution is a (slightly perturbed) rotation around y axis and the time averaged 
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Figure 6. Quantum fidelity in the perturbative regime i5 < for a = 30, J = 1500, 
and 7 = 7r/2 and 7r/6, calculated as a full trace Hilbert space average as a function of 
the scaled variable J{St)^. For 7 = tt/2 data ior S — 1 ■ 10^^ (solid curve) and 5 • 10~^ 
(dotted curve) are shown. For 7 = 7r/6, 6 = 3- lO"'^ (sohd) and 1 • 10"^ (dotted) 
are shown. Note that for both 7 the curves for both 6 practically overlap. The chain 
curves are theoretical predictions (p^) with classically computed ad- 



We will now use (approximate) analytical results for a — > to compare with numerics 
for a = 0.1. It should be noted that our leading order analytical approximations could 
easily be systematically improved using a classical perturbation theory (treating a as a 
perturbing parameter). However, since the agreement, as shown below, is almost perfect 
in all cases, we see no need for refinement at this level. 

First, we consider the full trace average, p = 1/J, and starting from the expression 
(p!5|), F{t) = (exp {itA6/h)), write the fidelity as a sum over all eigenvalues of Jy, namely 
(2m — 1)^, for m = 1, . . . , J/2 (in OE subspace). 



perturbation can be computed analytically (in the leading order in a) 




(53) 




J/2 



(54) 



m=l 



For large J we can replace the sum with an integral and get 



where erfi(2;) = Jq^ e ^ dt is a complex error function with a limit 

limj-^oo |erfi(ie*'^/^v^)l = 1 to which it approaches by oscillating around 1. We there- 
fore have an analytic expression for the fidelity (^Sj) in the case of an uniform average 




(55) 
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Figure 7. Fidelity in quasi-integrable regime for a = 0.1, 7 = 7r/2, S = 0.01, J = 100, 
and p — 1/J, in OE subspace. Solid curve gives the result of a numerical simulation. 
Isolated dots denote differences between numerical calculation and the analytic formula 
(HI) for a ||-Fnum.(i)P — |^anaii.(i)P|. The dotted line gives a predicted asymptotic 
decay cx t^^^^ , and the dashed/chain curves are the predicted fidelity decays at small 
times, namely the second order expansion = 1 — ( Ji(5)^/180, and 'improved' by 

the gaussian (p7[). 



over whole Hilbert space or, equivalently, over one random initial state. Its asymptotic 
decay is which agrees with the general semiclassical asymptotics p3| ) and we ex- 
pect initial quadratic decay (|T3|) for small times t < r^^ Decay rate t^c is determined 
by the time averaged correlation C (0) which can be calculated explicitly in the limit 
a — s> where the classical correlation function Cci(t) alternates for even/odd times as 
Cei(2t) = {z^{0)P{2t))/A = -1/90, Cd(2t + 1) = (52(0)z2(2t + l))/4 = 1/45, giving 

^^''"=° = K"^ + i5) = 1^' ^^^^ 



and the fidelity is expected to decay as (|T3|) with r^^ = vT80/(J5), for short times, 
t < Tne- The short-time formula (|13D and the full analytic expression (^5D are 
compared with the numerical simulation in fig. |^. The agreement is very good and, 
surprisingly enough, the gaussian approximation (|1^) for small times is observed to be 
valid considerably beyond the second order expansion (|1^). Quite interesting is the 
regime where the decay time t^c = VI8O/ {J5) for a "regular" dynamics will be smaller 
than a decay time rem = l/('^^^^o"ci) (|5TD for a "chaotic" dynamics. This will happen 



for 6 < 6r = l/(Jo"ci a/180) (|^). This border has the same scaling with J as 5mix 

Second, considering SU(2) coherent initial state (P^D, we could proceed along the 
same line, namely by an analytic calculation. But rather than that, we will illustrate 
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Figure 8. Fidelity for 7 = 7r/2, S = 0.025 and J = 100 on OF subspace. The dashed 
curve is a simulation for a = 30 (mixing regime, full trace average). The pluses are for 
a pure coherent initial state (see text for details) at a = 0.1 (non-mixing regime). The 
chain and solid curves are, respectively, the theoretical exponential (j^) and gaussian 
(|9|) decays. 



the usefulness of a semiclassical formula for F^2^{t) (^). This is a more general 
approach, as the explicit analytical calculation is usually not possible. Let us by -(9, if 
denote the spherical angular coordinates measured with respect to the y-axis. Then 
(/ = cos-i? = represent canonical action-angle coordinates for the integrable case 
a ^ 0. Furthermore, the coherent state (^) acquires a semiclassical gaussian form (p^) 
in the EBK basis \n), Jy\n) = n\n), namely 

|(n|^,V;)|^(xexp ^ . ), n = -. (57) 

The squeezing parameter A reads 

A = l/sin^t? = 1/(1 -y2). (58) 

In order to apply the general formula (^) we need to express the classical time average 
(^) in terms of a canonical action, a(J) = (1 — /^)/4, and evaluate the derivative, 
\da{l)/dl\'^ = giving \F^f{t)\ = exp {-6^ Jt^P {1 - P)/ 16). Rewriting this 

expression in terms of original spherical angles and (f measured with respect to z- 
axis, we obtain 

F;;f{t) = exp |^-^^{sin2 2^ sin^ ^ + sin^ ^ sin^ 2ip}^ . (59) 

Here the same interesting question, namely when do we have Tne-coh < ^em, results 
in a condition (|4^) 6 < 1/(4(TciV2J3) = 4/ V32crci In J (we did not write the factor 
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involving trigonometric functions (|59| ) which is of order 1). This condition generally 
can not be meet if we also require 6 > 6s (in order to see full exponential decay in 
the mixing regime), except if 32(Tciln J < 1. This is for instance satisfied for J = 100 
and 7 = 7r/2, for which = 0.00385. We checked this by a numerical simulation 
and the result for a coherent state centered at {'d,'^) = 7r(l/\/3, l/\/2), for which 
sin^ 219 sin^ + sin^ sin^ 2(^9 = 0.96, is shown in fig. |^. We can see that for t > 50 
the fidelity in a non-mixing regime (rne-coh = 23) is lower than the fidelity in a mixing 
regime (rem = 42). For larger times, t > t*, non- mixing decay F^°^ displays revivals of 
fidelity. 

Finally, we want to visualize the phenomenon of faster decay of fidelity in a regular, 
non- mixing regime a = 0.1 in comparison with the chaotic, mixing regime a = 30 by 
using the phase space representation of wave functions. A convenient and popular choice 
is a Husimi function H{ip, cos'i?) of a state {ip) defined as 

H{ip,cos^) = \{^,ip\^)\\ (60) 

We have chosen a pure random initial state {ip) and propagated it, in the first case 
with the propagator U for a = 0.1, and in the second case for a = 30, as well as with 
the perturbed propagators Us in both cases. The state and all the other parameters 
(e.g. perturbation strength S) were identical in both cases. Then we compared the 
differences between the Husimi functions AH{ip, cos{}) = H{ip, cos{}) — Hs{(p, cos-d) of 
the perturbed and the unperturbed time evolution. The results for a = 0.1 are shown 
in fig. ^ and for a = 30 in fig. |T0[ For chosen 6 = 10~^ the fidelity decay for the regular 
dynamics is much faster than for the chaotic one and this effect can also be observed 
in Husimi functions by comparing the right columns of both figures. But one should 
stress that by considering the difference of Husimi functions AH = H — Hs, we loose 
information about the relative phases of the perturbed and unperturbed wave functions 
in the coherent-state basis which are more important for the fidelity decay in the regular 
case (a = 0.1). However, results shown in figs. suggest that not only the quantum 
phases but also the amplitudes (in some classically motivated, e.g. coherent state basis) 
exhibit larger susceptibility to system perturbations for the regular (non-mixing) as 
compared to the chaotic (mixing) dynamics. 



3.2. Pair of coupled kicked tops 



As we have already remarked at the end of sect. |2.3| , for a one dimensional d = 1 
systems, the 'surprising' behavior Tne-coh < Tem is for coherent initial states possible 
only around the border (^I]) 6s (unless cTci is very small like in the example of fig. |^) 
where the exponential decay in the mixing regime goes over to a gaussian due to finite 
size Af. However, for more than one degree of freedom, d > 1, such behavior is generally 
possible well above the finite size - perturbative border 6s. In order to illustrate this 
phenomenon we will now briefly consider a numerical example of a pair of coupled kicked 



tops where d = 2. 
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Figure 9. Husimi function i/((y5, cos z9) for a regular evolution a = 0.1, 7 = 7r/2 
and three different times t — 0,500 and 1000 (left column). Gaussian random initial 
state is used and J — 200, S = 10^"*. In the right column we show pictures of 
absolute difference between Husimi functions of the unperturbed and perturbed time 
evolution, \H{(f,cos'd) — (95, cos i9) | . Fidelity at the times shown is |F(500)| ~ 0.73 
and |F(1000)| — 0.25. Top picture in the right column gives the classical phase space 
portrait. 



We consider coupled kicked tops with a unitary propagator (a simplified version of 
that of ref. H) 

U{e) = e~*f ■^i5'e~*^'^^^e~*'^'^^^'^^^'''^. (61) 

where J\ and J2 are two independent quantum angular momentum vectors. Perturbed 
propagator is obtained by perturbing the parameter e, so that Us = U{e + 6). The 
perturbation generator is therefore 

A = JizJ2z, (62) 

with h = 1/ J. A pair of coupled kicked tops possesses the same pair of discrete 
symmetries as the single top for 7 = n/2 0, namely and Ry, and in addition 
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Figure 10. The same as in fig. ^ except for a chaotic (mixing) dynamics, a — 30. 
Nmiierical values of the fidelity are now |F(500)| = 0.997 and |F(1000)| = 0.988. 



it is unvariant under the permutation of the identical tops. However in our simple- 
minded numerical experiment reported here we have used the propagator (|6T|) over the 
full (2J + 1)^ dimensional Hilbert space thus making the appropriate average over all 
symmetry classes. 

The classical limit is obtained by J ^ cxo and writing the classical angular 
momentum vectors in terms of two unit vectors on the sphere ri 2 = Ji,2/J- In 
component notation we get the following equations of motion 

2 = ^1,2 (63) 
y'i,2 = yi,2 cos(e2;2,i) + Xi,2 sin(e2;2,i) 
^1,2 = - ^1,2 cos(ex2,i) + 1/1,2 sin(e2;2,i) 

We have chosen two regimes, namely non-ergodic (KAM) regime for e = 1 where the 
vast majority of classical orbits are stable, and the mixing regime for e = 20 where no 
significant traces of stable classical orbits were found and where the correlation sum was 
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Figure 11. Fidelity for two coupled kicked tops, 5 = 8- 10^'* and J = 200. The 
upper curves are for e = 20 (mixing regime), solid curve for a coherent initial state 
( |65| ) and dashed curve for a random initial state, and the lower - dotted curve is for 
e = 1 (non- mixing regime) with a coherent initial state (|65|). The exponential and 
gaussian chain curves give, respectively, the expected theoretical decays (^Tj) and ( p7| , 
here decay time is determined by the best fit). 



to a very good accuracy given by the first term only 

1 1 /. 1 

2 

Our motivation here was to compare the non-ergodic and mixing fidehty decays for the 
coherent initial state which is here the dyadic product of SU(2) coherent states ( ^9|) 

1^9, <^)i2 = 1^92,^^2)2® |t9i,<^i)i. (65) 
In fig.|TT| we show the fidelity decay at J = 200 and 6 = 8-10"'^ in non-ergodic and mixing 
cases started from the coherent state (|65|) with = = 7r(l/-\/3, l/-\/2). 

We find excellent agreement between the theortical predictions (0) and (|27|) and the 
numerics. Note that we are here already in the regime 5 < 5c where the fidelity decay 
in ergodic-mixing regime is slower rem = 673 than in non-ergodic regime Tne-coh = 379. 
In the ergodic-mixing regime (e = 20) we show for comparison also the fidelity decay 
for a random initial state which is (due to ergodicity) almost identical to the case 
of coherent initial state. We note that overall fidelity decay here is similar as in a 
one-dimensional case, however, the scaling of various time and perturbation scales on 
h = 1/ J is different as discussed in subsect. |?^ . We want to conclude this section with an 
interesting experimental observation, namely that in the non-ergodic regime the revivals 
of fidelity (quantum reccurrences) beyond the time scale t^ are much less pronounced 
in 2d than in Id, e.g. compare figs.|^ and |ll]. 
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4. Classical fidelity 

In this section we wisli to contrast ttie surprising findings on the stability of quantum 
motion with an application of the general idea (section 2) to a fully classical concept, 
namely to (unitary) Perron-Frobenius propagator of classical phase space densities. Let 
us consider a classical dynamical system being given by an invertible measure (volume) 
preserving map <p : x ^ (pi^) on a D (= 2d) dimensional phase space Ai. No additional 
properties of the map will be assumed except of being at least piece-wise differentiable 
so that the stability matrix 

M(.) ^ (66) 

OX 

exists almost everywhere. Unitary Perron-Frobenius propagator Ud is defined on a phase 
space density from p G L^(A^) as 

iUdp)ix) = p{cl^^-'\x)). (68) 

The next step is to perturb the map in a most general way by composing it 

05 = o gs, (69) 

with a near identity map gsi^x) which is generated by an arbitrary (smooth) vector field 
a{x) thru the 5— flow 

= a(gs), with initial conditions go(x) = x. (70) 

do 

Vector field a{x) should be divergence free 

diva(a;) = (71) 

to keep the perturbed map (f)s volume preserving (or hamiltonian for symplectic maps). 
Now we can write the perturbed unitary Perron-Frobenius propagator 

{Ud,sp){x) = pig^f'\cl,^-'\x))) = (f/,iexp {-5 aix) ■ d/dx) p) {x) (72) 

in the standardized general form (|l|) with the self-adjoint perturbation operator A^i 

d 

Uc\,s = Uci exp(-«(5 Aci), Aci = -la- —. (73) 

We define the 'classical fidelity' as the inner product between the phase space 
densities propagated by two slightly different Perron-Frobenius propagators from the 
same fixed initial density p{x) 

Fd{t)= j d^x {Ui,^,p){x){Ul,p){x) = [ d^xp(0i-*)(=r))p(0(-*)(x)).(74) 
Jm Jm 

Note that -Fci(t) follows the corresponding quantum fidelity F[t) (2) for very short times 

t -C ts, since the latter can be alternatively written as F{t) = J d'^xW{x,t)Ws{x,t) 

in terms of Wigner functions W{x,t) and Ws{x,t) of the states U^lip) and U^\ip), 

respectively, corresponding to the Liouville densities t/^p and ^/cuP C^^) if p{x) = 
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W{x,t = 0). Fci(t) can be given a probabilistic interpretation, namely if p{x) is a 
characteristic function on a set i3 C A^, p{x) = 1 or 0, for x G or a:; ^ Ai, 
respectively, with the volume V = Jj^^x p{x) then V~^Fci(t) is a probability that a 
perturbed system with the initial condition chosen at random from the set B is found, 
after time t, in the image of the same set B propagated for time t with the unperturbed 
system. 

In order to proceed formally along eqs. we define and straightforwardly 

compute the 'classical Heisenberg operators' 

:= f/dMeit/*i = -zat{x) ■ ^, at{x) := M;\x)a{cf>^'\x)). (75) 

Finally, plugging this into 5— expansion (^, observing that the first order always vanishes 
since diva^ = 0, and simplifying the second order by integrating by parts, we obtain 



Fei(t) = 1-^1 I d^'x 

I Jm 



dp{x) 



t-i 



dx 



Y^m-,\x)a{ct>^-''\x)) 



2 



0{5^). (76) 



Note that now the RHS does not look like a correlation function at all! The reason 
for this lies in an essentially different form of the generator of perturbation ( [73D which 
in classical case involves phase space derivatives invoking the stability matrix, while in 
quantum mechanics it is a simple quantization of the observable a{x). 

Now, observe that the classical linear response formula (|76D gives precisely the same 
result as one intuitively expects based on the stability of individual orbits. Indeed, in 
the case of chaotic dynamics, the integrand (for sufficiently long times t') becomes 
dominated by the leading eigenvalue of the stability matrix M^, namely the largest 
Lyapunov exponent A, and (|76|) can be written as 

i^ci-cha(t) = 1 - consti (5^ exp(2At) + 0{5^) (77) 

where the unspecified constant depends on the initial distribution p{x). This means 
the classical fidelity for chaotic dynamic decays on logarithhmically short time-scale 
~ log(l/5)/A. On the other hand, in the case of regular dynamics, all eigenvalues of 
the stability matrices Mj have modulus 1 so the integrand remains bounded by some 
constant and the fidelity can be estimated as 

Fei_reg(t) = 1 - cousta 5H'^ + 0{5^) (78) 

Note that the fidelity decay for regular dynamics has the same 5 and t dependence as 
in the quantum chase (|13]) with the decay time-scale ~ 5"^. 



5. Summary and discussion 

In this paper we have presented a simple theory on intrinsic decoherence of deterministic 
(isolated) quantum dynamical systems due to small static perturbations of the evolution 
operator. The central object of study is the fidelity F{t) of quantum motion computed 
either with a single (pure) initial state or averaged over an ensemble of initial states 
described by the statistical density matrix. 
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The main result of the paper is a simple linear response (or Kubo-like) formula 
(^) which relates the fidelity decay to the total sum (or integral) of two point time 
autocorrelation function of the generator of perturbation. In the limit of infinite Hilbert 
space dimension we found exponential fidelity decay on a time-scale r^ra {h/6Y, for 
quantum ergodic and mixing systems, whereas for non-ergodic systems we have found 
much faster decay (in the sufficiently 'quantum' regime where 6 <^ h) on a time-scale 
(X h/6 for random initial states or Tne-coh oc Ti^^"^ /5 for coherent initial states (minimal 
uncertainty wavepackets) where the fidelity F{t) is given by a Fourier transformation 
of the local density of states of the time averaged perturbation operator. A special 
emphasis was given to the semiclassical theory of fidelity of small but finite values of ^, 
where different regimes and the corresponding time and perturbation scales are carefully 
discussed, and where fidelity decay may asymptotically (as ^ — >• 0) be evaluated in terms 
of classical quantities only. Interestingly, finite size fiuctuations of fidelity (for very long 
times at a finite Hilbert space dimension) have been shown to be given by the inverse 
participation ratio of the eigenstates of the perturbed evolution operator in the eigenbasis 
of the unperturbed propagator. The surprising aspects of our relations are mainly due 
to non-interchangability of the limits 5-^0 and — > as the relevant decay time- 
scales are only functions of the ratio 5/h. Therefore, a different and intuitively expected 
behaviour, namely faster fidelity decay for mixing than regular dynamics, is obtained 
in the 'classical' regime where % <^ 5 (or making the limit — > prior to 5 ^ 0). 
A similar, reassuring result has been found by applying our formalism to inspect the 
analogous classical fidelity for the unitary Perron-Frobenius evolution of volume (area) 
preserving maps: there the classical fidelity for regular dynamics has been found to 
decay on a time-scale oc while for a chaotic dynamics, the fidelity decay is governed 
by a maximal Lyapunov exponent A on a short time-scale oc ln(l/5)/A. 

We note that our findings should be of primary importance for understanding the 
problem of stability of quantum computation |Tl[] where different regimes and time-scales 
^ should now be explained in terms of intrinsic dynamics of a particular quantum 



algorithm. Actually, our relation between fidelity and correlation decay may prove 
very useful in design of any new technology which is aimed at manipulating individual 
quantum states. On the other hand, our results may also shed new light on the relation 
between decoherence and dynamics In particular, due to quite special role of 



coherent initial states on the short Ehrenfest time-scale one may expect the results to 
be quite different for a random initial state and/or for longer time-scales (see |P7|| ). 
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